Tin telluride: a weakly co-elastic metal 
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We report resonant ultrasound spectroscopy (RUS), dilatometry/magnetostriction, magnetotrans- 
port, magnetization, specific-heat, and ^^®Sn Mossbauer spectroscopy measurements on SnTe and 
Sno.ggsCro.oosTe. Hall measurements at T = 77 K indicate that our Bridgman-grown single crystals 
have ap-type carrier concentration of 3.4 x 10*® cm~^ and that our Cr-doped crystals have an n-type 
concentration of 5.8 x 10^^ cm~^. Although our SnTe crystals are diamagnetic over the temperature 
range 2K < T < HOOK, the Cr-doped crystals are room temperature ferromagnets with a Curie 
temperature of 294 K. For each sample type, three-terminal capacitive dilatometry measurements 
detect a subtle 0.5 micron distortion at ~ 85 K. Whereas our RUS measurements on SnTe show 
elastic hardening near the structural transition, pointing to co-elastic behavior, similar measure- 
ments on Sno.ggsCro.oosTe show a pronounced softening, pointing to ferroelastic behavior. Effective 
Debye temperature, 6d, values of SnTe obtained from **®Sn Mossbauer studies show a hardening of 
phonons in the range 60-115K {6d = 162K) as compared with the 100-300K range (Od ~ 150K). 
In addition, a precursor softening extending over approximately 100 K anticipates this collapse at 
the critical temperature, and quantitative analysis over three decades of its reduced modulus finds 
AC44/C44 = A\{T — ro)/To|~'^ with K = 0.50 ± 0.02, a value indicating a three-dimensional soften- 
ing of phonon branches at a temperature To ~ 75 K, considerably below Tc. We suggest that the 
differences in these two types of elastic behaviors lie in the absence of elastic domain wall motion 
in the one case and their nucleation in the other. 



I. INTRODUCTION 

For around 50 years, the IV- VI compounds have pro- 
vided a basis for semimetals and small-gap semicon- 
ductors with high-temperature phases in various crys- 
tal structures transforming into other structural phases 
on lowering the temperatur^^. The low-temperature 
phases themselves have often been objects of consider- 
able interest. On structurally transforming, SnTe and 
GeTe, for example, become polar, with some authors 
invoking ferroelectricity. In practice the physical prop- 
erties of members of this family of compounds depend 
on the value of their excess carrier concentratiorPEl, al- 
though an experimental link to defects is also possible. 
In fact, several ternary compounds, such as Pbi-^^Suj^Te 
and Gei_j;Sna;Te, have also received considerable study 
because changing the carrier concentration was seen as 
way not only to tailor the size of the band gaps but also 
to switch the parit y of th e bands across the gap. A re- 
cent interest revivaP*- ** * *^ * in these materials is occurring 
in part because of a perceived opportunity to create novel 
tunable multifunctional devices on the nano-scale. Just 
as recently, we started a study of SnTe to exploit another 
property of this and several other IV- VI compounds; 
namely, upon doping with a magnetic atom, such as Cr or 



Mn, the resulting compound becomes ferromagnetic^^. 
Thus, magnetically-doped SnTe has the potential of be- 
ing multiferroic with three features - ferroelasticity, fer- 
roelectricity, and ferromagnetism - in a considerably 
simpler crystal structure than the perovskite-structured 
manganite multiferroics that typically are studied for just 
two features - ferroelectricity and ferromagnetism. It was 
envisaged that the combination between semi-conducting 
and ferroic properties is desirable to investigate the pos- 
sibility of a new class of ferroic materials. 

In this paper, we report the presence of ferroelasticity 
in Cr-doped SnTe below the phase transition point near 
lOOK. Ferroelasticity implies that a spontaneous strain 
arises below a critical temperature and that the switching 
of the variants must be allowed. In contrast, SnTe is co- 
elastic, i.e. the appearance of spontaneous strain occurs 
without variant switching, where the elastic strain energy 
is insufficient to form ferroelastic domains. 

More specifically, the high-temperature phase of SnTe 
transforms from cubic (rocksalt) to rhombohedral around 
100 K, the precise value depending on the excess car- 
rier concentration. Electronic structure theory points to 
unsaturated bonding dominated by the p-orbitals of Sn 
and Te with a stoichiometric compound being insulating 
with an approximate 2 eV gap across the fee Brillouin 
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zone except a round the L-point where the gap is about 
0.3 eV^EHMll Across the L-point, the band states have 
opposite parity and couple to a transverse optical phonon 
in the [111] direction, a Raman active mode that drives a 
second-order displacive transition in this direction. The 
measured size of the displacement is approximately 50 
pm. Accompanying it is a shear strain of approximately 
0.122° that decreases the rhombohedral angle. 

Prima facie, the transition in SnTe seems to fit the 
scenario discussed by Anderson and BlounllA^, whereby 
a second-order transition with strain must be accompa- 
nied by some other internal symmetry change. One pos- 
sible change is the spontaneous development of a symme- 
try represented by a polar (vector) order parameter. In 
the present case, this polar order parameter is the static 
displacement associated with the optical phonon. Such 
a polar order parameter would produce ferroelectricity 
if there were a spontaneous electrical polarization that 
can be reversed by the application of an electric field. 
The best evidence for electrical polarization in SnTe is 
its anomalous transverse effective charge inferred from 
Raman scattering^. For the observed sizes of the dis- 
placement, the effective charge implied by the measured 
dipole moments is around 4 to 8 times the unit charge 
available from fully ionized Sn and Te. Another cited 
feature is the large static and optical frequency dielec- 
tric constants of approximately 1000 and 50, respectively, 
near the structural transition. 

In this paper we report a new suite of measurements 
on SnTe that produced unexpected results. In practice 
SnTe is always found to have a substantial carrier density, 
probably from defect-induced carriers, and we confirmed 
this in our samples. In consequence an internal electri- 
cal field cannot be maintained and ferroelectricity as a 
macroscopic phenomenon is not possible. We did ob- 
serve weak signatures of a transition, presumably struc- 
tural, around 100 K; however, instead of ferroelasticity 
we found co-elasticity. Across the transition, ferroelas- 
tic and coelastic crystals exhibit only small changes in 
physical properties, such as elastic stiffness and thermal 
expansion, whereas ferroelastics typically exhibit large 
changes. The co-elastic changes we found for SnTe are 
likely among the smallest yet observed. 

We also performed measurements on Cr-doped SnTe 
and found several unexpected results. Specifically, we 
found a 0.5 % Cr doping produced a room temperature 
ferromagnet emerging from diamagnetism in the high- 
temperature phase. More surprisingly, the dilute level of 
doping changed the material from co-elastic to ferroelas- 
tic. 

The purpose of this paper is to report these findings. 
We remark that our samples, as other samples of SnTe re- 
ported in the literature, are off-stoichiometry, and hence 
metallic. Elsewhere,!^ we have reported angular-resolved 
photoemission (ARPES) measurements that show over a 
wide range of temperatures a Fermi surface whose topol- 
ogy changes markedly with electron binding energy. A 
hole doping of about 5 x lO^^cm"^ shifts the fermi en- 



ergy into the uppermost valence band by several tenths 
of an eV, consistent with our magnetotransport results. 
In this paper we focus on other measurements that sup- 
port the presence of the transition, the metallic behavior, 
and the weak elastic and thermal anomalies. We remark 
that the strong carrier-concentration dependence of the 
structural transition in SnTe is an effect normally not 
seen in other structurally transforming materials, such as 
shape-memory alloys where the martensitic phase tran- 
sition is structurally identical to a ferroelastic transition 
but where the equivalent Landau free energy change is 
very larg^i^. 

It is convenient to begin from a theoretical perspective 
because we will find it straightforward to interpret some 
of the results using a Ginzburg-Landau theory. In Sec- 
tion II, we outline a basic dual order-parameter Landau 
free-energy functional for studying the phase transition. 
In Section III, we discuss how we grew our SnTe crys- 
tals and how we performed our resonant ultrasound spec- 
troscopy (RUS), dilatometery, magnetostriction, magne- 
totransport, specific heat, and ^^^Sn Mossbauer spec- 
troscopy measurements. In Section IV, we report the 
results of our measurements. Finally, in Section V, we 
give summary remarks. 



II. THEORY 

We will discuss the structural phase transitions in pure 
SnTe and Cr-doped SnTe in terms of a Landau- Ginzburg 
free energjff^HlSl^ As standard, we construct this en- 
ergy by starting from the symmetry of the two phases. 
The high temperature phase has the rock salt structure 
(space group Fm3m, site symmetry m3m, Wykoff posi- 
tions 000 and (^ ^ ^)), and the low temperature phase has 
the rhombohedral structure (space group R3m, site sym- 
metry m-3m, Wykoff positions (xxx) with x = 0.004 and 
i - 0.004, rhombohedral angle 59.878 at zero K'j^. By 
group theory arguments, the phase transition occurs by a 
two step symmetry reductiorP2l. The first is the reduction 
of the cubic phase symmetry to that of a rhombohedral 
phase (R-3m) via an elastic strain along the cube diag- 
onal. The irreducible representation of this elastic order 
parameter is T2^^, and this resulting phase has an inver- 
sion center. The allowed strain components are xy, xz, 
and yz. The second symmetry reduction removes the in- 
version center by a relative displacement of the two sub- 
lattices whose irreducible representation is Ti^^^. The 
basis functions of Ti„ are x, y, and z. 

The primary order parameter describes the polar dis- 
placement between the two sublattices of approximately 
50 pm. (In SnTe, this represents a one part in 10^ effect.) 
This displacement vector is parallel to (111). The sec- 
ondary order parameter is the strain \exy \ = \exz\ = \s-yz\, 
leading to an elongation of the unit cell along any of the 
polar axes. This shear corresponds to an rhombohedral 
shear angle of 0.122°. The conjugated elastic modulus 
is C44, which is expected to be explicitly temperature 
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dependent in the vicinity of the transition, whereas all 
other elastic moduli remain relatively unaffected by the 
phase transitiorfiSl. Both order parameters act at the ori- 
gin of the Brillouin zone and maintain the topology of 
the unit cell. 

The Landau approach to this transformation leads to 
a mean-field approximation where all secondary interac- 
tions are treate d as perturbations of the primary order 
parametej ^^ l ^^ l. Band structure calculation^ report a 
non-convex energy well that is deep along the polar de- 
formation but is shallow along the spontaneous strain, 
a situation typical for multiferroic materials where elec- 
trical and elastic degrees of freedom mb^. We will fol- 
low the treatment of ferroelectric materials and treat the 
polar order parameter as the driving (primary) order pa- 
rameter and add the elastic relaxation in the appropriate 
symmetry-adapted perturbation. The Gibbs free energy 
of the zone center cubic-to-rhombohedral transition is 

G{Pi, Sjk) — L{Pi) + /cl(ejfe) -I- /coupling (A J Gjk) 

+ /Gi„zburg(VP,) , (1) 

where the first term on the right is the non-convex Lan- 
dau energy, the second term is the elastic energy, the 
third term is the coupling energy between the polar or- 
der parameter Pi and spontaneous strain e^fc components, 
and the last term is the Ginzburg energy which is rele- 
vant for the assessment of the domain boundaries. In 
the following paragraphs, we state expressions for only 
the first three of these contributions, thereby choosing to 
describe the physics in a single domain. 

The general form of the Landau fre e energ y L{Pi) for 
the rhombohedral (mSm) symmetry iJEESESl 

LiP) \a{Pl + Pi + Pi) + \B'{Pt + Pl + Pi) 

2 D 

+ ^-C"{Pl + Pi + Pi)^ 

+ \c"'{Pl + P| + Pi){Pt + Pi + Pi) . (2) 

Within one domain, for example, the (111) domain, one 
component of Pi is non-zero while all others are zero. 
Thus the mono-domain Landau free energy has the re- 
duced form 

L(P) = laP' + -BP^ + icp6 , (3) 
2 4 5 

with coefficients B = B{B' , B") and C = C(C', C", C"). 
This is the simplest Landau free energy for a continuous 
phase transition when P > or P = and C > 0. The 
case P = is called tricritical as it separates a second 
order (P > 0) and first order (P < 0) transformation. 

The standard assumption in a Landau theory of phase 
transitions is a = A\T — Tc\. Very close to Tc, this as- 
sumption imposes a mean-field critical exponent on the 
thermodynamics, which is generally incorrect. Very far 



from Tc, as T — 7> 0, the temperature dependences of ther- 
modynamic quantities become inconsistent with the re- 
quirement that the entropy must approach zero. In the 
interpretation of our results, we will assume that the 
structural transition is driven by a single anharmonic 
soft mode. With this single-mode assumption, the low- 
temperature form 

a = AOs [coth(e,/r) - coth(e,/T,)] , (4) 

exactly follows fr om a self-consistent phonon approxima- 
tion to this mod^i^lSU. In various contexts, this form has 
provided informative fits to low-temperature properties 
derived from a Landau free energy. 

In Eq. Q, the parameter Qg is called the saturation 
temperature. SnTe has a Raman active soft mode in 
the rhombohedral phas^2Zl_ in this case, the saturation 
temperature is an estimate of its frequency. With these 
soft optical-mode frequencies-^-. We estimate the sat- 
uration temperature of SnTe to be between 40 K and 
60 K, which is consistent with the value obtained from 
our thermal expansion data (see Fig. |3|. The phononic 
saturation temperature based on a Debye temperature 
would be appropriate for a proper ferroelastic system. 

For the [111] domain the spontaneous strain is symme- 
try constrained with 64 = e^y + &xz + ^yz ■ The elastic en- 
ergy is for a single domain crystal /ci(e4) = iC44e|, and 
the coupling energy is /coupling (Pj, eij) = \P'^e^. If we 
were to relax the free energy with respect to this strain, 
that is, solve dG/de^ — and compress ([2]), the resulting 
expression would have the same functional form as Eq. [Sl 
but with P replaced by a renormalized parameter B*^. 
Thus, as typical when a secondary order parameter is 
present, the nature of the transition remains described 
by the Landau free-energy form of the primary cyrder pa- 
rameter but with a shifted critical temperature^-. 

The corollary to the effect of strain coupling in shift- 
ing the transition, is that a spontaneous strain appears 
in response to the order parameter. Because this is a 
macroscopic effect on the crystal, it is often much easier 
to observe than the order parameter itself. 



III. EXPERIMENTAL 

We prepared large single crystals (1x1x2 cm^) of SnTe 
either by fusing the elements in a Bridgman furnace (tem- 
perature gradient). We took advantage of the high vapor 
pressure of Te and lowered the growth temperature in 
the Bridgman to 850° C. Powder diffraction and back- 
reflection Laue images showed that the crystals grew 
along the (321) crystallographic directions, and rocking 
curves measured by neutron diffraction indicated a small 
mosaic of 0.5°. 

We determined the elastic constants of disk-shaped 
SnTe single crystals with resonant ultrasound spec- 
troscopy (RUS)r^. In this method we measure the me- 
chanical resonances of a sample held lightly between two 
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piezoelectric transducers, and deduce the complete elas- 
tic stiffness tensor from the spectrum of resonant frequen- 
cies by using a computer fitting algorithm. In our RUS 
apparatus, the sample and the transducers are mounted 
inside a home-built measurement cell, which in turn fits 
inside a Quantum Design Physical Properties Measure- 
ment System (PPMS). We excited the sample using a 
PZT transducer (Valpey Fisher) and a Stanford Research 
Systems (SRS) DS345 function generator and then mea- 
sured the sample response using a PZT transducer and 
an SRS SR844 lock-in amplifier. A Cernox thermome- 
ter (LakeShore Cryotronics) , located immediately adja- 
cent to the sample, gave us the sample temperature. 
Custom- written Lab View software controled the function 
generator and the lock-in and also instructed the PPMS 
firmware to measure the temperature and magnetic field. 
For our RUS measurements, frequency scans (200-1000 
kHz) occurred at approximately 60 temperatures in the 
range 300-20 K. Before starting each scan, we allowed the 
sample temperature to stabilize at the pre-programmed 
value, and we did all measurements in an approximately 
1 torr He exchange gas to ensure rapid thermalization 
between the sample, transducers, and surroundings. 

To check for evidence of SnTe piezoelectricity, we per- 
formed a modified RUS experiment. We did these ex- 
periments exactly as discussed above, except instead of 
using a PZT transducer as the pick-up transducer, we 
used a SnTe specimen. Specifically, with silver epoxy we 
attached the center conductor of a mini-coax cable to one 
side of a SnTe specimen and the shield of the coax to the 
other side. Then, we connected the other end of the coax 
to the input of the lock-in amplifier, exactly as we did 
for the PZT pick-up transducer. Finally, we placed the 
SnTe crystal between the two PZT transducers. One of 
the PZT transducers excites the sample, while the other 
holds it in place. Performed this way, our RUS uses the 
SnTe samples as pick-up transducers to self-detect their 
own resonances. 

We also made thermal expansivity/striction, specific 
heat, magnetic susceptibility, and Hall measurements 
using the PPMS. We measured the specific heat by a 
thermal-relaxation techniqu^^, and we used a vibrating 
sample magnetometer (VSM) with a frequency of 20 Hz 
and a large bore coilset to measure the magnetic suscep- 
tibility in fields up to 9 T. To measure the linear thermal 
expansion of SnTe and Sno.ggsCro.oosTe single crystals, 
we used a capacitive dilatometei^^ that fits inside the 
PPMS. All expansion measurements were done in He ex- 
change gas as the sample slowly warmed from 4 K to 
300 K at 0.2 K min~^. The capacitance was measured 
using a three terminal capacitance bridge with a reso- 
lution of 10~^ pF, which corresponds to a sensitivity of 
approximately 0.003 A. 

The material also lends itself to analysis with the 
Mossbauer effect, because of the ^^^Sn in SnTe. 
Mossbauer studies of the temperature and pres- 
sure dependencies can provide site specific informa- 
tion on the isomer shift (valence), the magnetiza- 



tion, the electric field gradient, lattice dynamics, phase 
transitions, local structure, and mixed composition. 
Mossbauer spectroscopy was carried out with ^^^SnTe 
and ^^^Sno.ggsCro.oosTe samples in the 65-300K range 
using a Ca^^^^SnOs source. A top-loading cryostat was 
used in which both source and absorbers were kept at the 
same temperature. A separate measurement in which a 
standard ^^^Sn absorber was kept at ambient temper- 
ature. The source temperature was varied in the 50- 
300K range. This allowed to extract the effective De- 
bye temperature, Od = 300i^, of Ca^^^SnOg needed 
for the analysis of SnTe lattice dynamics. Measure- 
ments with the ferromagnetic Sng.ggsCro.oosTe sample 
showed no measurable magnetic hyperfine interaction at 
the -'^^^Sn nucleus to the lowest temperature. 



IV. RESULTS AND DISCUSSION 
A. Transport, specific heat, and magnetization 

Hall effect measurements indicated that our pure SnTe 
samples had a p-type carrier concentration of 3.4x10^^ 
cm~'^ and our Sng.ggsCro.oosTe crystals had an n-type 
carrier concentration 5.8x10^^ cm~^. Based on tran- 
sition temperature vs. carrier concentration data in 
Ref. 1281 we expected the structural transition in pure 
SnTe to occur at aproximately 100 K. The apparent n- 
type doping in the Cr- doped material is anomalous, and 
we will return to it later. 

To locate T^, more precisely, we first performed resis- 
tivity and specific heat measurements. Our resistivity 
measurement vs. temperature data showed two unam- 
biguous breaks in slope, one at approximately 82 K and 
the other at 98 K. While the latter was near the expected 
critical temperature, the former was unexpected and we 
have no explanation for it. In general, on heating and 
cooling over this temperature range we were unable to 
retrace the resistivity curve. The measurements also de- 
pended on how and where we attached our leads. As we 
cooled to 6 K, the resistivity dropped rapidly to zero, in 
a manner consistent with the onset of superconductivity. 
Indeed, there have been previous reports of low tempera- 
ture, c arrier concentration dependent superconductivity 
in SnTePaUtall. 

Below 6K our specific heat measurements were fitted 
by Cp = jT + l3T^, with values of 7=0.32mJK-2mor ^ 
and the inferred Debye temperature, &d = 206K, values 
consistent with previous report^^. Using our measure- 
ments of 7 and the magnetic susceptibility, we computed 
the Wilson ratio and found it to be approximately 2, 
the free electron value. This result prompted two other 
types of measurements. One type was ARPES. As dis- 
cussed at length elsewher^^l^ both above 98 K and below 
82 K these measurements clearly indicated the presence 
of a Fermi surface that is consistent with ab initio band 
structure calculations. 

For our second type of measurement, we first tried 
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ferroelectric switching of SeTe (and Sno.ggsCro.oosTe). 
We found that any electric field applied to the sam- 
ples was short-circuited by their high electric conductivi- 
ties. Thus, while this short-circuiting excludes the direct 
measurement of ferroelectricity as a ferroic property, we 
tested whether the weaker condition of piezoelectricity 
could still apply. Our experimental results clearly show 
that this is not the case: On structurally transforming 
the elastic stresses do not generate dipole moments. A 
literature search failed to find previous reports of piezo- 
electric measurements. 

From these measurements we concluded that our non- 
stoichiometric SnTe is a metal, but we still needed to lo- 
cate the phase transition. In the vicinity of the expected 
critical temperature we detected no jump in the specific 
heat. This result is consistent with one previous reporP^ 
but inconsistent with anotheif^. Our ARPES measure- 
ments however did indicate that a transition likely occurs 
but the observed changes were too small to infer the sym- 
metry of the low temperature phase. X-ray diffraction 
measurements indicated a likely transition around 98 K 
but again out experimental resolution was inadequate to 
determine the symmetry of the low temperature phase. 
From these additional measurements, we assumed that 
a transition indeed occurs around 98 K, that it is struc- 
tural, and that the symmetry of the low temperature 
phase is rhombohedral. 

Sno.ggsCro.oosTe was found to show a ferromagnetic 
transition at around room temperature, with a low tem- 
perature moment consistent with full magnetization of 
the Cr. It is surprising to find such a high transition 
temperature in a lightly doped material; we remark that 
the compound Cr2Te3 is a ferromagnet with a Curie tem- 
perature of 325 K. Although we have not observed phase 
segregation on even microscopic scales, we believe that 
this is an important hint which will help to resolve re- 
sults that follow. 



B. Elastic properties 

We turn now to the clastic properties, which show clear 
evidence of a structural transition. 

Figure [3] shows the temperature dependence of the 
spontaneous strain in the R3m phase of a randomly 
oriented Sng.ggsCro.oosTe single crystal. The sponta- 
neous strain is defined as [AL/L] spontaneous = AL/L — 
[AL/L]b, where AL/L is the measured thermal strain 
and [AL/L]b is the baseline thermal strain if no phase 
transition occurred. We determined [AL/L]b by extrap- 
olating the expansion for the high-temperature Fm3m 
phase into the temperature range of the low-temperature 
R3m phase. To do this we assumed that the expansion 
of the high-temperature phase arises mainly from the an- 
harmonicity of the soft phonon mode that is associated 
with the phase transition. We therefore modeled the 
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FIG. 1: Temperature evolution of magnetization of 
Sno.ggsCro.oosTe. In this plot, the Curie- Weiss region is only 
in the temperature range for T > 300K. The parameters from 
a fit to the data above this region gives the Weiss constant 
Tw ~ 290K, indicating ferromagnetic ordering. Above the 
ordering region, the effective number of Bohr magnetons is 
p = 3.3. For paramagnetic Cr^"'" and Cr'^"'', we have p = 4.9 
and for p = 3.8, respectively. Thus p — 3.3 is a reasonable 
value, and indicates the Cr is in the -1-3 valence state. 



high-temperature expansion using an Einstein oscillator. 



[AL/L]b = D + F coth {9 m/T) , 



(5) 



where 9m is the characteristic temperature of the soft 
mode and D and F are constants. We found the value 
of 9m using an iterative procedure in which we guessed 
a value, fitted Eq. ([5| to the AL/L data for the high- 
temperature phase (in the range 90 < T < 145 K) to de- 
termine D and F, and then computed [AL/L] spontaneous- 
Using different values of 9m, we repeated this procedure 
until (i[AL/L]spontanoous/rfr — ^ at K, as required 
by the third law of thermodynamics. In this way, we 
found 9m = 45 K. This value is in good agreement with 
the soft-mode frequency measured using Raman spec- 
troscopjB^. As discussed in the theory section, this spon- 
taneous strain is associated with the elastic shear modu- 
lus C44. 

The remaining contribution is the spontaneous strain, 
which is shown in Fig. [3] For comparison, in Fig. [2] 
we report results of Mossbauer measurements. In the 
harmonic approximation using the Debye approximation, 
the recoilless fraction / can be written as: 



/ = exp 




(6) 



where is the ^Sn recoil energy (2.57 

meV=30.8K/A:B). In the high temperature limit, 
we have 



, , 6ErT 



T>l9n 



(7) 
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FIG. 2: Temperature dependence of the recoilless fraction, /, 
for the temperature regimes: (a) 60-115 K and (b) 90-300K, 
respectively. The derived 9d values., of 162K and 150K, re- 
spectively, indicates a hardening of the phonon spectrum in 
the low-temperature phase. 



The temperature dependent area A{T), of the absorption 
spectrum is proportional to the recoilless fraction of the 
sample /a(T): 



A{t) = const ,fa{T) , 



(8) 



where the constant includes all experimental parameters 
independent of T. The Debye temperature, ^d, can be 
extracted from the T-linear dependence of ln(/Q). In 
Fig. [2] we show the T-dependence of A for two temper- 
ature regimes, namely, the 60-115 K and the 90-300K 
regimes, respectively. Correspondingly, we obtain the 6d 
values as 162K and 150K, respectively, which indicates a 
hardening of the phonon spectrum in the 60-115 K versus 
the 100-300 K range. 

These results are consistent with the Landau the- 
ory discussed in the previous section. First, we de- 
rived a thermodynamic potential G{P) by relaxing the 
strain under the condition dG{P, Cij) /deij = where P 
is determined by the condition of thermal equilibrium 
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FIG. 3: (Color online) Temperature evolution of the sponta- 
neous strain, [AL/LJapontanooua, iu the low-temperature R3m 
phase. The analysis used to determine the excess strain is 
discussed in the text. The curve through the data is the best 
fit of Eq. ( 10 1. Parameters deduced from this fit are given in 
the text. 



dG{P)/dP = and e^- scales proportional to P^. We 
then used our experimentally measured thermal expan- 
sion for Sno.ggsCro.oosTe to determine the parameters of 
the Landau potential. 

For the present purposes, we simplify to treat the po- 
larisation and spontaneous strain as scalars, so the free 
energy becomes 



G{P,e) = ^aP'- 



1 



1 



1 



/3P^ + --fP^+SeP^ + -Ke^-ea (9) 
6 2 



Here we choose the model form of Eq[4] for the quadratic 
coefficient a. For convenience, we scale the polarisation 
to be dimensionless, by setting 7 — a{T = 0). (In the 
case of a tricritical theory (/3 — 0) will fix P{T = 0, e = 
0) = 1). e is the dimensionless strain, a the applied 
stress, and K the relevant elastic modulus. In this case all 
the material constants a, /3, 7, 6 and K have dimensions 
of energy per unit volume. 

The ground state is determined from the simultaneous 
minisation of Eq. ^ with respect to strain e and polar- 
isation P. In the absence of applied stress, this leads to 
a strain e = AL/L cx P^ which becomes 



Aa6s [coth(6'^/r) 



where b 



(^)( 



) and 



AS 



coth(^^,/T,)]) , 

(10) 

. The curve in 



Fig.[3]shows the best fit of (10 1 to the spontaneous strain 
data for SnTe and Sno.ggsCro.oosTe. From this fit we de- 
duced a = 1.09 X 10-® K-i and 6 = 6.3 x 10"*, 6^ = 
72.1 K, and = 82.7 K. Notice that b'^/Aa « 9 K, indi- 
cating that the phase transition is close to tricriticaPE^. 
This analysis allows us to extract the coupling constant 
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to strain, namely 6/K = 5.8 x 10 ^, and the renormalised 
quartic coefficient P' h = (/3 - 25'^IK)h = 1.08. 

We checked the vahdity of these parameters by calcu- 
lating the temperature evolution oi P — Po^/eA where 
Po is a proportionality factor. The calculated values of 
P were compared with the experimental one^^ and the 
good agreement affirms that our normalized Landau po- 
tential is a good approximation for the excess Gibbs free 
energy of the transition. 

We next estimated an upper bound for the absolute 
values of the Landau parameters. A direct determina- 
tion would require experimental data for the excess spe- 
cific heat of the phase transition in order to calibrate the 
specific heat jump ACp = T (fG{P)/dT^\T^. Hatta and 
KobayashP3 measured the excess specific heat at temper- 
atures near the transition point, and they found a jump 
at Tc of 0.45 J mor^K"^ We were unable to reproduce 
their data. Our own resolution was below 0.2Jmol-iK-i 
but we still could not detect a reproducible signal mark- 
ing the transition in pure SnTe. If we take the value 
of 0.1Jmol~"'^K~^ as an estimate for the excess value at 
Tc, we can estimate an upper bound for A, specifically, 
A = 0.13 J mor^K-^ 

It has been known for some time that both increasing 
carrier concentration and increased pressure suppress the 
transition, though reliable pressure data has only been 
obtained on a Pbi-xSuxTe alloy. With the estimated 
parameters we can now predict the pressure-dependence 
of the critical temperature. Under an externally imposed 
stress, there is a Hooke's law strain of e = a/K, which 
contributes a extra positive quadratic term in the polari- 
sation, and therefore shifts the transition temperature to 
the temperature where a + 26e = a{T) + 2Sa/ K vanishes. 
The critical pressure as a function of temperature is then 



Pc{T) 



A0.. 



2{S/K) 



coth(^) -coth(^) 



(11) 



With the above parameters, the transition is predicted 
to vanish beyond a pressure of 1.9kbar, shown in Fig- 
ur^ We are not aware of measurements on pure SnTe, 
but Sugai et have reported results on the alloy 

Pbo.56Sno.44Te, which has a lower critical temperature 
and pressure. This data is also shown in Figur^ fit to 
the same functional form as ( 11 1, with quite clearly sim- 



ilar parameters to those we extracted for the higher- T^, 
material. 

With a further assumption we can interpret the car- 
rier concentration dependence from previous measure- 
ments, also shown in Fig. |4] We note that we can 
analyze the carrier concentration dependence n of the 
phase transition with the same perturbation approach. 
The symmetry-allowed lowest-order term for the car- 
rier concentration dependence of G{P) is kP^ n so 
for a continuous transition the condition for Tc becomes 
a* = a + kn = or Tc(n) = e,5/(coth"^(coth(es/Tc(n = 
0)) - kn)). In Fig.|4) we compare the calculated depen- 
dence of the transition temperature on the carrier con- 
centration, with k = 0.37 X 10^° cm~'^, with experimental 
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FIG. 4: (Color online) Dependence of the transition temper- 
ature on pressure and on the carrier concentration as derived 
from Landau theory including the quantum saturation (data 
points from Sugai et al.l^Sl). Shown also is the predicted pres- 
sure dependence of Tc from the Landau parameters estab- 
lished in this paper. 



datcP. 

A change in the elastic properties is indicated near the 
transition, and we now turn to direct measurement of 
the shear modulus from resonant ultrasound. Figure |5] 
shows the temperature dependence of [/i//i(300 K)] , 
where /i is the lowest mechanical resonant frequency of 
the sample measured using RUS, and /i(300 K) is its 
value at 300 K. Data are shown for pure SnTe, measured 
at zero applied magnetic field, and for Sno.ggsCro.oosTe, 
measured at and 250 Oe. For both pure and Cr-doped 
SnTe, the change in /i depends only on the value of C44 
(for example, dfi/dC44 = 1 whereas dfi/dC — and 
dfi/dB — 0). Thus Fig. [5] shows the temperature depen- 
dence of the normalized shear modulus C44/C44 (300 K). 
Our results are consistent with Beattie's observation^ 
that pure SnTe has no elastic softening at the transition, 
but rather a small stiffening in the R3m phase below ap- 
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FIG. 6: (Color online) Log-log plot of the precursor soften- 
ing versus the reduced temperature T — Tq. The theoretical 
softening is AC44 = A{T - To)"" with k = 0.5. 



below Tc the internal friction grows in the Cr-compound 
and not in the pure material. This makes it clear that 
this is not the effect of critical fluctuations, but likely 
a precursor efl-ect due to nucleation and motion of do- 
main walls. In Fig. [6] we show that the elastic precursor 
softening AC44/C44 of Sno.ggsCro.oosTe fits to a power 
law in reduced temperature T — Tq, where Tg is a tem- 
perature some 10 K below the critical temperature for 
the phase transition. The precursor softening is the dif- 
ference between the measured temperature dependence 
of (744/(744(300 K) and the baseline T-dependence if no 
phase transition occurred. For the baseline we used the 
temperature evolution of (744/(744(800 K) in the high- 
temperature phase of pure SnTe, extrapolating these data 
linearly for temperatures below 110 K so as to exclude 
the phase transition. 

The elastic precursor effect in Cr-doped SeTe is much 
larger than that usually seen in ferroelastic material^^. 
The temperature dependence of the precursor softening 



given by theory is 



AC44/C44CX [(r-To)/To] 



(12) 



proximately 100 K. In contrast, we find that for Cr-doped 
SnTe there is a dramatic softening of C44 associated with 
the transition. The softening is accompanied by a large 
increase in the damping of the mode, shown in the lower 
panel of the figure, absent for the pure compound. 

We note that well above the transition temperature 
(for example, T > 175 K), C44 shows the same temper- 
ature scaling in both the pure and the Cr-doped mate- 
rials, yet softening in the Cr-doped material occurs well 
above Tc. 

Note that the softening observed in the Cr-doped ma- 
terial extends from well above the phase transition tem- 
perature (Tc w 85K) to well below it. Notice further that 



where the value of k depends of the dimensionality of 
the phonon softening. For the high-temperature phase 
of Sno.ggsCrg.ogsTe, with cubic space group Pm3m, we 
would expect a value oi k ~ 0.5, consistent with 3- 
dimensional softening of the phonons. This is not a criti- 
cal phenomenon, but represents a three-dimensional soft- 
ening of the acoustically coupled soft mode. 



V. CONCLUDING REMARKS 

The transport and (lack of piezoelectricity) measure- 
ments reported here and the recently reported ARPES 
measurements are sufficient to establish that SnTe is a 
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metal and not a semiconductor. Hence, it is not a ferro- 
electric as commonly assumed. Without doubt the ob- 
served metallic behavior, instead of the semiconducting 
behavior predicted by band structure calculation, is a 
consequence of the non-stoichiometry. We see this result 
as significant because we are unaware of any published 
claim of having prepared this material stoichiometrically. 

We note a number of our measurements produced in- 
conclusive results about the existence, symmetry, and lo- 
cation of the phase transition. We believe this is partially 
a consequence of the subtlety of the transition. Although 
we have no direct evidence, we also believe some mea- 
surements, in particular, those of the low temperature 
elastic constants, were influenced by inhomogeneities in 
the material. Likely, SnTe is not mono-domained. In 
this regard, we remark that our ultrasound measurements 
confirmed Beattie's observation^^ that low-temperature 
SnTe does not soften elastically but instead hardens. The 
absence of softening indicates that mobile twin bound- 
aries do not nucleate at the transition point. Such be- 
havior is extremely unusual, having apparently so far 
been observed only in martensitic Cu74.08Al23.13Be2.79 
where the immobility of the twin boundaries is caused by 
strong dislocation pinnin^^. All other martensitic trans- 
formations studied to date show strong elastic softening 
near T^, similar to what we observed in Suq ggsCrg oosTe 
(Fig.[5]). 

Several physical scenarios can explain the observed 
hardening in SnTe. One is that the twin boundaries 
are heavily pinned by dislocations as in the case of 
Cu74.o8Al23.i3Be2.7g^. In the present case, however, this 
scenario is unlikely because by symmetry the phase tran- 
sition in SeTe does not generate dislocations, and other 
pinning mechanisms for domain walls, such as Peierls 
pinning, are too weak to prevent the movement of fer- 
roelastic domain boundaries. A second more probable 
scenario is the absence of domain boundaries. This sce- 



nario implies what we measure at the transition is the 
intrinsic hardening of the shear modulus. In contrast to 
SnTe, we note that Cr-doped SnTe near exhibits the 
strong elastic softening which we suggest is related to 
the nucleation of ferroelastic domain boundaries and to 

their mobility under the externally applied stress of the 
ultrasouncP^l, 

We noted earlier that our Cr-doped samples show bulk 
ferromagnetism at a temperature close to that expected 
for the compound Cr2Te3. We have seen no microstruc- 
tural evidence for phase segregation, but it is likely that 
a nano-scale second phase plays an important role in 
nucleating domain boundaries and enhancing their mo- 
bility. Cr2Te3 is a good metal, and metal inclusions 
in the semimetallic SnTe would also produce anomalous 
magnetotransport. For practical purposes, however the 
Cr-SnTe samples show homogeneous ferromagnetism and 
ferroelasticity. Coupling between the magnetic and struc- 
tural order is nonetheless weak, and probably indirect. 

In closing, we comment that while future experiments 
may clarify the material science of the elastic properties, 
it is interesting to propose the more fundamental ques- 
tion, "If stoichiometric SnTe existed, would it be ferro- 
electric at low temperatures?" If so, then some amount of 
carrier doping, perhaps, an infinitesimal amount, would 
generate a very unusual metal-insulator transition, that 
is, one accompanied by the creation or destruction of 
electrical polarization. 
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